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Pteridine reductase (PTR1) is a potential target for drug development against 
parasitic Trypanosoma and Leishmania species, protozoa that are responsible 
for a range of serious diseases found in tropical and subtropical parts of the 
world. As part of a structure-based approach to inhibitor development, 
specifically targeting Leishmania species, well ordered crystals of L. donovani 
PTR1 were sought to support the characterization of complexes formed with 
inhibitors. An efficient system for recombinant protein production was prepared 
and the enzyme was purified and crystallized in an orthorhombic form with 
ammonium sulfate as the precipitant. Diffraction data were measured to 2.5 A 
resolution and the structure was solved by molecular replacement. However, 
a sulfate occupies a phosphate-binding site used by NADPH and occludes 
cofactor binding. The nicotinamide moiety is a critical component of the active 
site and without it this part of the structure is disordered. The crystal form 
obtained under these conditions is therefore unsuitable for the characterization 
of inhibitor complexes. 
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1. Introduction 

Leishmania are protozoan parasites (order Trypanosomatida, class 
Kinetoplastida) that cause a range of diseases and present a serious 
health risk to millions of people worldwide (Desjeux, 2004; 
Reithinger et ai, 2007). The incidence of infection is primarily in 
tropical and subtropical regions of the world (Herwaldt, 1999). The 
different species of Leishmania are responsible for distinctive 
conditions (Reithinger et ai, 2007). For example, L. donovani causes 
visceral leishmaniasis, a potentially fatal disease, while infection with 
L. major leads to cutaneous leishmaniasis. Several compounds are 
available to treat these infections, but increasing levels of drug 
resistance combined with the high cost and toxicity of antileishmanial 
drugs compromises the control of the diseases (Croft et ai, 2006; 
Maltezou, 2010). These observations explain in part why the World 
Health Organization has identified the leishmaniases as neglected 
diseases and is urgently seeking novel therapeutic approaches (World 
Health Organization, 2007). 

Our aim is to identify and characterize drug targets in these 
parasites and to apply structure-based approaches to develop potent 
inhibitors that possess the right chemical properties to underpin 
early-stage drug discovery (Hunter, 2009). A promising target with 
respect to infection with Leishmania sp. is the NADPH-dependent 
short-chain dehydrogenase/reductase pteridine reductase (PTR1; EC 
1.5.1.33). This enzyme is unique to trypanosomatid parasites, where it 
supports the provision of reduced biopterins that are necessary for 
metacyclogenesis (Cunningham et al., 2001) and which are implicated 
in resistance to reactive oxygen and nitrogen species in Leishmania 
(Moreira et al, 2009). PTR1 catalyzes the reduction of biopterin to 
7,8-dihydrobiopterin as well as its subsequent reduction to 5,6,7,8- 
tetrahydrobiopterin. Additionally, the enzyme catalyzes similar 
reactions in the salvage of unconjugated folates in Leishmania 
(Nare, Hardy et al., 1997; Nare, Luba et al., 1997). As Leishmania are 
auxotrophic for pteridines (folates and pterins) and are required to 
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obtain these nutrients from their environment to maintain growth, 
disrupting this salvage process represents a potential therapeutic 
strategy. 

We have previously studied the structure-mechanism-activity 
relationships for the enzymes from L. major (LraPTRl) and Tryp- 
anosoma brucei (TZ?PTR1) and determined the structures of a series 
of inhibitor complexes (Gourley et al, 2001; Dawson et aL, 2006; 
Cavazzuti et aL, 2008; Mpamhanga et aL, 2009), whilst the structure of 
T. cruzi PTR2 has been reported by others (Schormann et aL, 2005). 
We identified sequence and structural differences between LraPTRl 
and LZ?PTR1 that explain why some inhibitors display a significant 
level of selectivity for one orthologue over the other (Gibson et aL, 
2009; Tulloch et aL, 2010). Although we were able to routinely 
generate crystals of 77?PTR1 that diffracted to between 2.0 and 1.0 A 
resolution (Dawson et aL, 2010), studies with LraPTRl have been 
hampered by poor crystal quality and a lack of reproducibility. One 
crystal form of LraPTRl diffracted to beyond 2.0 A resolution but 
could only be obtained in the presence of NADPH and methotrexate 
(Gourley et aL, 2001); when other ligands were present different 
crystal forms were obtained. The size of the asymmetric unit is 
increased from two to either four or eight subunits and the crystals 
are often mechanically twinned and diffract to lower resolution, 
with the diffraction pattern being anisotropic and highly mosaic 
(McLuskey et aL, 2004; Schiittelkopf et aL, 2005). An alternative 
source of Leishmania PTR1 was therefore sought for our investiga- 
tions. Studies with L. tarentolae PTR1 resulted in a 2.8 A resolution 
structure of the complex with NADPH, but despite its presence in the 
crystallization mixture the tight-binding ligand methotrexate was not 
observed in the electron-density maps (Zhao et aL, 2003). This was 
not considered to be an improvement on the results that we had 
previously obtained, so we elected to initiate crystallographic studies 
of L. donovani PTR1 (LdPTRl), this also being the enzyme from the 
pathogen that causes the most serious form of leishmaniasis. There is 
a high level of conservation (91% sequence identity) between the 
L. major and L. donovani enzymes and homology modelling of the 
latter has suggested a close structural relationship in and around the 
active site (Kaur et aL, 2010). 



2. Methods 

2.1. Expression and purification 

The gene encoding LdPTRl was cloned into the expression vector 
pET15b (Novagen) modified to encode a tobacco etch virus (TEV) 
protease-cleavable N-terminal hexahistidine tag. This plasmid was 
heat-shock transformed into Escherichia coli BL21 (DE3) GOLD 
cells (Stratagene) and selected on Luria-Bertani (LB) agar plates 
containing 50 mg l -1 carbenicillin. Cells were then cultured in LB 
medium containing the same antibiotic at 310 K with shaking. Gene 
expression was induced with isopropyl /3-D-l-thiogalactopyranoside 
at a final concentration of 1 mM when the cells reached the mid-log 
phase of growth (optical density at 600 nm of 0.6-0.8). After incu- 
bation at room temperature for 16 h, the cells were harvested by 
centrifugation (4000g, 277 K, 30 min), washed with fresh medium and 
the centrifugation was repeated. Cell pellets were frozen at 253 K 
until required. 

For purification, the frozen cell pellet from 1 1 culture was thawed 
on ice and resuspended in lysis buffer (50 mM Tris-HCl pH 7.7, 
200 mM KC1, 20 mM imidazole) supplemented with an EDTA-free 
protease-inhibitor cocktail tablet (Roche) and 100 |ig DNAse I 
(Sigma- Aldrich). Cells were lysed using a French cell press and the 
lysate was clarified by centrifugation (37 OOOg, 277 K, 30 min). The 



supernatant was filtered (0.2 um) and applied onto a 5 ml HisTrap 
metal-chelating column (GE Healthcare) preloaded with Ni 2+ and 
equilibrated with lysis buffer. A linear imidazole gradient was applied 
and the protein eluted at a concentration of approximately 160 mM 
imidazole. Fractions containing LdPTRl were pooled and the histi- 
dine tag was cleaved by incubation with TEV protease at 303 K for 
3 h. Imidazole was removed from the buffer by dialysis and cleavage 
continued at 277 K overnight before a second Ni 2+ -affinity step was 
performed, isolating the pure and cleaved LdPTRl. The protein was 
concentrated to approximately 10 mg ml -1 using a centrifugal filter 
unit with a molecular-weight cutoff of 3500 Da (Millipore) and the 
buffer was exchanged to 20 mM Tris-HCl pH 7.7. The high level of 
purity (>95%) and the molecular weight of the protein were con- 
firmed by SDS-PAGE and matrix-assisted laser desorption/ionization 
time-of-flight mass spectrometry (data not shown). The yield of 
purified protein was relatively low at approximately 3 mg per litre of 
bacterial culture and we note that most of the material produced was 
insoluble. 

2.2. Crystallization and data collection 

Sitting-drop vapour diffusion was used to test a range of 
commercially available screens. Crystals were only obtained in the 
presence of both cofactor and an inhibitor and these were optimized 
in hanging drops. A solution of 5 mg ml" 1 LdPTRl, 1 mM NADP + , 
20 mM dithiothreitol and 1 mM methotrexate was incubated on ice 
for 1 h before crystallization. The protein solution (1 ul) was mixed 
with reservoir solution (1 ul) and stored at room temperature. Small 
orthorhombic plates with approximate dimensions of 0.1 x 0.1 x 
<0.05 mm formed within days above reservoirs containing 0.1 M MES 
pH 6.5, 10%(v/v) dioxane and 1.6 M ammonium sulfate (Fig. 1). 

Crystals were placed in a nylon loop and flash-cooled to 100 K in a 
stream of nitrogen gas after first being cryoprotected by passing them 
through a solution of 40% PEG 400. Initial in-house X-ray experi- 
ments using a Rigaku MicroMax-007 rotating-anode X-ray generator 
and an R-AXIS IV ++ image-plate detector produced only smeared 
low-resolution diffraction images. The quality of the diffraction was 
improved by an annealing step in which the cryostream was manually 
diverted and reintroduced after 5-10 s (data not shown). Diffraction 
data were then collected on beamline 104 at the Diamond Light 
Source synchrotron. 




Figure 1 

A crystal of LdPTRl. The approximate dimensions of this sample were 0.1 x 0.1 x 
0.03 mm. 
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2.3. Structure solution and refinement 

The data were processed and scaled using MOSFLM (Leslie, 2006) 
and SCALA (Evans, 2006; Collaborative Computational Project, 
Number 4, 1994), respectively. The space group is C222 x and the unit- 
cell parameters are a = 107.51, b = 126.44, c = 87.51 A. The Matthews 
coefficient (Matthews, 1968) of 2.5 A 3 Da -1 , which corresponds to 




Figure 2 

(a) Ribbon diagram of the LdPTRl tetramer. Subunits A and B constitute the 
asymmetric unit, (b) Ribbon diagram of an LdPTRl subunit. Seven red /3-strands 
(numbered) are sandwiched between the blue a-helices; loop regions are coloured 
yellow. A sulfate is depicted as sticks in the active site with S coloured orange and O 
red; the N- and C-termini of the protein are labelled. Blue, red and black asterisks 
mark residues either side of the missing @3-a3, @4-a4 and @6-a6 loops, 
respectively. 



Table 1 

Data-collection and refinement statistics. 
Values in parentheses are for the highest resolution shell. 
Space group 

Unit-cell parameters (A) 
Resolution range (A) 
Wavelength (A) 
No. of measurements 
No. of unique reflections 
Multiplicity 
Completeness (%) 
Mean Ila(I) 
Wilson B factor (A 2 ) 

^merget (%) 
^workt (%) 
tffree§ (%) 

R.m.s.d. bonds (A) 
R.m.s.d. bond angles (°) 
Ramachandran analysis 

Favoured (%) 

Allowed (%) 

Outliers (%) 
Protein residues (total) 
Atoms (total) 
Overall B factor (A 2 ) 
Additional groups 

Waters 

No- 
Average B factor (A 2 ) 
Sulfates 

No- 
Average B factor (A 2 ) 



C222 l 

a = 107.5, b = 126.4, c = 87.5 

29.9-2.5 (2.6-2.5) 

0.973 

144524 (21138) 
21004 (3022) 
6.9 (7.0) 
99.9 (100.0) 
11.9 (3.8) 
46.9 

9.6 (42.4) 
22.7 (28.0) 
28.5 (33.0) 
0.019 
1.802 

95.2 

4.6 

0.2 

432 

3201 

43.5 



24 
38.0 



2 

53.4 



t Emerge = EhkiZi IW - (Met)) I/£mz £/ I t (hkt), where Iihkl) is the intensity of 
the ith measurement of reflection hkl and (I(hkl)} is the mean value of Ii(hkl) for all i 
measurements. $ # work = J^hki \ \ F obs I - l^cak 1 1 /Em; \ F obs I . where ^obs is the observed 
structure-factor amplitude and F calc is the structure-factor amplitude calculated from the 
model. § R tree is the same as R wor ^ except calculated with a subset (5%) of data that 
were excluded from refinement calculations. 



approximately 50% bulk solvent, suggested that two subunits occupy 
the asymmetric unit. 

The LdPTRl structure was solved by molecular replacement using 
Phaser (McCoy et al., 2007). The search model was a single LraPTRl 
subunit (PDB code le7w; Gourley et al, 2001) in which divergent 
residues were truncated to and all ligands were removed. Two 
subunits were positioned and rigid-body refinement and all subse- 
quent refinements were carried out using REFMAC5 (Murshudov et 
al, 1997). The graphics program Coot (Emsley & Cowtan, 2004) was 
used to inspect difference and electron-density maps for model 
fitting, solvent and ligand searching. Sulfate and water molecules 
were added after the protein model was completed and further 
rounds of refinement were performed. Noncrystallographic symmetry 
restraints were not employed during the analysis. Crystallographic 
statistics are shown in Table 1. Figures were prepared with PyMOL 
(DeLano, 2002). 



3. Results and discussion 

3.1. Overall structure 

PTR1 is a tetrameric enzyme and this crystal form of LdPTRl has 
two subunits (labelled A and B) in the asymmetric unit; a 2 1 screw 
axis parallel to c generates the tetramer (Fig. 2a). Each monomer is 
formed by a seven-stranded central /3-sheet flanked by a set of three 
a-helices on either side (Fig. 2b): a Rossmann-fold repeat (Gourley et 
al, 2001). The structures of subunits A and B are highly conserved, 
with an r.m.s.d. of 0.51 A when 212 C* atoms are matched; therefore, 
unless otherwise stated all descriptions refer to subunit A of LdPTRl. 
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3.2. A disordered active site 

Sequence comparisons between LdPTRl and LraPTRl (data not 
shown) indicate that the residues involved in the construction of 
the active site and that are essential for catalysis in binding cofactor, 
substrates, products and inhibitors are strictly conserved and are also 
highly conserved in TZ?PTR1 (Dawson et al, 2006). However, the 
residues near the catalytic site were poorly ordered in contrast to the 
same regions of LraPTRl and 7WTR1, with no electron density 



Sulfate 




Figure 3 

The F Q — F c difference density OMIT map for sulfate bound in the LdPTRl 
NADP + -binding site (green mesh contoured at 3.5a). The orange dotted lines 
represent potential hydrogen bonds formed between the ion and the enzyme. These 
interactions are in the range 2.7-3.1 A. 



Methotrexate 




Figure 4 

A C a trace of one LdPTRl subunit (cyan) overlaid with an LraPTRl subunit (PDB 
code le7w; Gourley et al, 2001; black). The sulfate ion bound to LdPTRl is shown 
as yellow sticks, while the NADPH and methotrexate binding to LraPTRl are 
depicted as red and orange sticks, respectively. The @6-a6 loop of the LraPTRl 
model is marked. This loop is absent from LdPTRl. The @5-a5 substrate-binding 
loop adopts different positions in the two structures. 



corresponding to NADP + or methotrexate, ligands that were present 
in the crystallization conditions in a fivefold molar excess. The K { for 
methotrexate inhibition of LraPTRl is 39 ± 19 nM and that with 
respect to 7WTR1 is 152 ± 16 nM (Dawson et al., 2006). 

Ammonium sulfate was the precipitant for crystal growth and a 
sulfate ion binds in a polar cavity formed by the /31-cd turn and the 
loop between /32 and al, accepting hydrogen bonds donated from 
three main-chain amides (His38, Arg39 and Ser40) and Ser40 OG 
(Fig. 3). This polar cavity is the binding site for the adenine 
2'-phosphate group of the cofactor (Gourley et al., 2001; Dawson et 
al., 2006). 

The core structure of the subunit is preserved between the 
LdPTRl and LraPTRl structures (Fig. 4) and an overlay of 202 C 
atoms common to both subunits resulted in an r.m.s.d. of 1.3 A 
(calculated using Coot, Emsley & Cowtan, 2004). Owing to the 
absence of well defined electron density, the LdPTRl model contains 
fewer residues compared with those of LraPTRl or LZ?PTR1. Missing 
segments include residues 70-80 (the loop linking /33 and a3), 112- 
132 (the fi4-a4 loop) and 227-254 (the fi6-a6 loop). The fi3-a3 loop 
is also poorly ordered in LraPTRl (Schiittelkopf et al., 2005). 
However, the /34-a4 loop is well ordered in LraPTRl but is missing in 
LdPTRl. This means that in LdPTRl, Phell3, which is a critical 
residue, is disordered. The phenylalanine, together with the nicotin- 
amide, forms 7T-stacking interactions that stabilize ligand binding 
in the catalytic site (Gourley et al., 2001; Dawson et al., 2006). The 
largest break in the electron density involves residues 227-254, which 
form what is termed the substrate-binding loop linking /36 to a6 
(Tulloch et al., 2010). Residues in this loop bind parts of folate sub- 
strates, products and some inhibitors. The absence of the /34-a4 loop 
vacates an area in the active-site cleft, allowing the /35-a?5 loop of 
LdPTRl to adopt a different position to fill the gap. Here, the C a 
atoms of several residues relocate by between 10 and 16 A. 

In LdPTRl Argl7 is disordered. This residue is important in 
binding the cofactor pyrophosphate (McLuskey et al., 2004), which in 
turn interacts with and positions the nicotinamide. Aspl81, Tyrl94 
and Lysl98 are the key catalytic residues. Aspl81 is located within 
the link between /35 and a5, but this loop has been built into weak 
electron density relative to the structure as a whole. The orientation 
of Tyrl94 is similar to that observed in LraPTRl. However, only the 
C a atoms of Lysl98 and its closest neighbours agree reasonably well 
with the identical LraPTRl residues, while the side chain extends into 
a position that is inappropriate to form the stabilizing interaction 
formed with the nicotinamide ribose in LraPTRl (Gourley et al., 
2001; data not shown). 

PTR1 displays a sequential ordered mechanism, with first cofactor 
binding and then substrate; following reduction the product leaves, 
followed by oxidized cofactor (Luba et al., 1998; Gourley et al., 2001). 
Substrate or inhibitors can only bind after the binary protein- 
cofactor complex has formed. Since sulfate binding blocks the 
2'-phosphate-binding site, NADP + is not present and this leads to 
disorder in the substrate-binding part of the active site. 



4. Conclusions 

The structure of apo LdPTRl has been solved to 2.5 A resolution. 
Attempts to improve the LdPTRl crystal quality by attempting to 
crystallize the apo form and the binary complex with cofactor and by 
introducing other ligands in combination with oxidized and reduced 
cofactor failed. Even though we see no evidence for these ligands 
in the electron-density maps, the presence of both methotrexate and 
NADP + was essential to obtain this crystal form. Despite their 
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presence in the crystallization mixture, there was no electron density 
corresponding to these molecules. Instead, the high ammonium 
sulfate concentration in the crystallization conditions resulted in the 
replacement of an NADP + phosphate by sulfate ions. The formation 
of the catalytic site of PTR1 is dependent on the presence of nico- 
tinamide and in the absence of NADPH( + ) this part of the active site 
is disordered. We conclude that owing to the disorder this crystal 
form is not suitable for characterizing the interactions of LdPTRl 
with inhibitors. 

This research was funded by The Wellcome Trust (WT082596 and 
WT083481). We thank the Diamond Light Source for synchrotron 
beam time and excellent staff support. 
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